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Abstract 

The  nonlinear  development  of  the  most  unstable  Gortler  vortex  mode  in  boundary 
layer  flows  over  curved  walls  is  investigated.  The  most  unstable  Gortler  mode  is  confined 
to  a  viscous  wall  layer  of  thickness  and  has  spanwise  wavelength  it 

is,  of  course,  most  relevant  to  flow  situations  where  the  Gortler  number  G  ^  1.  The 
nonlinear  equations  governing  the  evolution  of  this  mode  over  an  0(G“®/®)  streamwise 
lengthscale  are  derived  and  are  found  to  be  of  a  fully  nonparallel  nature.  The  solution  of 
these  equations  is  achieved  by  making  use  of  the  numerical  scheme  used  by  Hall  (1988)  for 
the  numerical  solution  of  the  nonlinear  Gortler  equations  valid  for  0(1)  Gortler  numbers. 

Thus,  the  spanwise  dependence  of  the  flow  is  described  by  a  Fourier  expansion  whereas  the 
streamwise  and  normal  variations  of  the  flow  are  dealt  with  by  employing  a  suitable  finite 
difference  discretization  of  the  governing  equations.  Our  calculations  demonstrate  that, 
given  a  suitable  initial  disturbance,  after  a  brief  interval  of  decay,  the  energy  in  all  the 
higher  harmonics  grows  until  a  singularity  is  encountered  at  some  downstream  position. 

The  structure  of  the  flow  field  as  this  singularity  is  approached  suggests  that  the  singularity 
is  responsible  for  the  vortices,  which  are  initially  confined  to  the  thin  viscous  wall  layer, 
moving  away  from  the  wall  and  into  the  core  of  the  boundary  layer. 
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1.  Introduction 

Our  concern  is  with  the  effect  of  nonlinearity  on  the  growth  of  the  most  unstable  Gortler 
mode  in  the  high  Reynolds  number  flow  over  a  wall  of  variable  curvature.  The  most 
unstable  Gortler  mode  was  recently  identified  by  Denier,  Hall  and  Seddougui  (1991)  who, 
in  the  context  of  the  receptivity  problem  for  Gortler  vortices,  re-examined  the  linear 
Gortler  equations  in  the  limit  of  large  Gortler  number.  This  analysis,  while  considering 
both  the  inviscid  Gortler  vortex  modes  of  0(1)  wavenumber  and  the  right-hand  branch 
modes  of  wavenumber  (G  being  the  Gortler  number)  previously  considered  by 

Hall  (1982a),  was  able  to  demonstrate  the  existence  of  a  most  unstable  Gortler  mode. 
This  mode  was  found  to  lie  in  the  wavenumber  regime  0(G^/®)  with  a  spatial  growth 
rate  0(G^^^)  and  was  found  to  be  confined  to  a  thin,  0(G-i/5),  viscous  layer  located 
at  the  material  boundary.  This  should  be  compared  with  the  inviscid  and  right-hand 
branch  Gortler  modes  which  have  growth  rates  0(G^/^).  Note  that  it  is  readily  shown 
that  for  the  temporal  evolution  of  Gortler  vortices  the  growth  rates  are  of  0(G^/^)  over 
the  whole  wavenumber  spectrum;  in  fact  the  temporal  evolution  of  the  Gortler  vortex 
modes  is  of  no  practical  interest  since  in  situations  where  Gortler  vortices  are  found  to 
occur  experimentally  their  initial  development  is  apparently  independent  of  time. 
Previous  investigations  into  the  nonlinear  evolution  of  Gortler  vortices  are  limited  to  the 
work  of  Hall  (1988),  in  which  the  fully  nonlinear  Gortler  governing  equations  were  solved 
numerically  for  both  0(1)  wavenumber  and  Gortler  number,  and  that  of  Hall  and  Lakin 
(1988),  in  which  a  fully  nonlinear  mean  flow-first  harmonic  theory  is  developed  for  the 
right-hand  branch  Gortler  modes.  (See  also  Aihara  (1976)  who  derived  a  nonlinear  differ¬ 
ential  equation  to  determine  the  evolution  of  Gortler  vortices;  this  calculation,  however, 
ignores  nonparallel  effects  present  due  to  boundary  layer  growth  and  employs  other  ap¬ 
proximations  which  cannot  be  justified).  The  results  of  Hall  (1988),  which  are  applicable 
to  0(1)  wavenumbers  and  Gortler  numbers  for  which  a  full  numerical  treatment  of  the 
nonlinear  equations  governing  the  vortex  evolution  is  required,  demonstrated  that  as  the 
vortex  evolved  downstream  the  energy  of  the  flow  becomes  concentrated  in  the  mean  flow 
correction  and  fundamental  harmonic.  These  rc.-uiis,  coupled  with  the  fact  that  the  effect 
of  nonlinearity  in  this  problem  is  to  stabilize  the  flow,  suggest  that  the  far  downstream 
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structure  of  the  velocity  field  is  governed  by  a  first  harmonic/mean  flow  theory.  Moti¬ 
vated  by  these  results,  together  with  the  earlier  weakly  nonlinear  theory  of  Hall  (1982b), 
Hall  and  Lakin  (1988)  subsequently  developed  a  first  harmonic/mean  flow  theory  for  short 
wavelength,  large  amplitude  Gortler  vortices.  (We  note  that  the  small  wavelength  approx¬ 
imation  inherent  in  this  theory  is  not  physically  unrealistic  since  for  a  growing  boundary 
layer  the  effective  wavenumber  of  the  Gortler  vortices  increases  in  the  downstream  direc¬ 
tion).  This  theory  then  demonstrates  that  the  mean  flow  adjusts  due  to  the  presence  of 
the  vortex  state  so  as  to  render  the  large  amplitude  vortices  neutrally  stable.  The  mean 
flow  is  then  driven  by  the  vortex  velocity  field  and  its  form  is  completely  altered  from  that 
which  exists  in  the  absence  of  the  vortex  motion.  Subsequently  Hall  and  Seddougui  (1990) 
showed  that  the  large  amplitude  vortex  states  found  by  Hall  and  Lakin  (1988)  are  unstable 
to  wavy  vortex  modes  trapped  in  the  shear  layers  bounding  the  region  of  vortex  activity. 
The  above  results,  however,  are  not  directly  applicable  to  situations  where  the  Gortler 
number  is  large  and  the  wavenumber  is  not  close  to  its  neutral  value.  Here  we  would  expect, 
on  the  basis  of  the  linear  receptivity  theory  of  Denier  et  al  (1991),  the  most  unstable,  linear, 
Gortler  mode  to  be  excited  by,  for  example,  small  localized  surface  imperfections.  One  such 
area  where  this  scenario  would  be  applicable  is  the  flow  over  turbine  blades  where,  due 
to  the  large  Reynolds  number  encountered  in  such  situations  together  with  an  appreciable 
curvature,  we  would  anticipate  large  Gortler  numbers  and  thus  the  possible  excitation  of 
the  most  unstable  Gortler  mode  through  the  presence  of  small  surface  defects.  For  these 
reasons  a  study  of  the  nonlinear  evolution  of  the  most  unstable  Gortler  modes  is  both 
warranted  and  is  of  real  practical  importance. 

The  outline  of  the  rest  of  the  paper  is  as  follows.  In  §2  we  derive  the  nonlinear  equations 
governing  the  downstream  evolution  of  the  most  unstable  Gortler  vortex  mode.  A  brief 
discussion  of  the  numerical  tchcme  employed  to  solve  these  equations  is  given.  In  §3  we 
present  the  results  of  our  numerical  calculations  along  with  a  discussion  of  these  results. 
Finally  in  §4  we  draw  some  conclusions. 

2.  Formulation  of  the  governing  equations 

Consider  the  flow  of  a  viscous  incompressible  fluid,  of  density  p  and  kinematic  viscosity  u, 
over  a  wall  of  variable  curvature  a~^  Here  a  and  L  are  the  typical  length  scales 
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associated  with  the  radius  of  curvature  of  the  wall  and  the  downstream  development  of  the 
flow,  respectively.  Denoting  by  Uoo  the  free  stream  speed,  sufficiently  far  from  the  wall, 
we  define  the  usual  Reynolds  number  Re  by 

Re  =  (2.1) 

u 

and  restrict  our  attention  to  the  limit  Re  — ►  oo  with  the  Gortler  number,  G,  defined  by 

OT 

G  =  — Rei/2,  (2.2) 

a 

held  fixed.  Let  (x^y^z)  denote  the  usual  Cartesian  co-ordinates,  nondimensionalized  with 
respect  to  L,LRe~^^^ ,LRe~^^^  respectively.  The  corresponding  velocity  field  is  taken  to 
be 

u  =  Uooiu  +  u,  Re~^/^{v  +  u),  Rt~^^^w),  (2.3) 

where,  we  will  assume  that,  (u{x,y),v{x,y))  corresponds  to  the  Blasius  boundary  layer 
and  (u,  V,  w)  and  the  corresponding  pressure  perturbation  p  are  functions  of  (x,  y,  z).  By 
substituting  the  above  relations  into  the  Navier-Stokes  and  continuity  equations  we  find 
the  governing  equations  for  u,v,w  and  p,  correct  to  G(Re“^/^),  to  be  given  by 

Wr  +  +  ICz  =  0, 

Uyy  +Uzz  —  VUy  =  UUx  +  UxU  +  VUy  +  Ql  , 

(2.4) 

Vyy  +  Vzz  -  GXUU  -  Py  =  UVx  +  +  VVy  +  VyV  +  Q2, 

Wyy  +  Wzz  —  Pz  =  UWx  +  VWy  +  Q^. 

Here  the  functions  Q\,Q2,  Qa  are  nonlinear  and  axe  given  by 

Ql  =  UUx  +  VUy  +  WUz, 

Q2  =  UVx  +  VVy  +  wvz  +  ^Gxu^,  (2.5) 

Qa  =  UWx  +  VWy  +  WWz- 

By  setting  the  nonlinear  terms  Qi,Q2,Qa  equal  to  zero  in  (2.4)  we  recover  the  linearized 
equations  of  Hall  (1933).  In  the  limit  of  large  spanwise  wavenumber,  d/dz  >  1,  Hall 
(1982a)  demonstrated  that  in  this  regime  the  right-hand  branch  of  the  neutral  curve 
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scales  as  G  =  0(k‘^),  where  A:  1  is  the  spanwise  wavenumber,  and  that  the  evolution  of 

the  Gortler  vortices  in  this  regime  can  be  readily  understood  without  recourse  to  any  type 
of  ad  hoc  assumptions  concerning  the,  so  called,  nonparallel  terms  which  are  inherently 
present  in  the  system  due  to  the  underlying  boundary  layer  growth.  A  subsequent  re¬ 
analysis  of  the  linearized  form  of  the  system  (2.4)  by  Denier,  Hall  and  Seddougui  (1991) 
demonstrated  that  in  the  large  Gortler  number  limit,  G  >•  1,  there  is  a  Gortler  vortex 
mode  of  maximum  growth  rate  G(G^/®).  It  is  found  to  occur  in  the  wavenumber  regime 
k  =  G(G^/®)  and  is  confined  to  a  thin  viscous  layer  situated  at  the  wall.  This 

should  be  compared  with  both  the  right-hand  branch  calculation  of  Hall  (1982a),  who 
found  that  the  growth  rate  for  this  mode  is  of  order  G^^^  with  wavenumber  0(G'/^)  and 
is  confined  to  a  thin  internal  viscous  layer  and  the  inviscid  Gortler  vortex  mode  considered 
by  Denier  et  al  (1991),  which  was  found  to  have  growth  rate  0(G^/^)  with  wavenumber 
G(l);  see  figure  1  for  a  schematic  representation  of  the  wavenumber /growth  rate  curve. 
For  0(1)  values  of  the  Gortler  number  the  fully  nonlinear  system  (2.4)  has  been  solved 
numerically  by  Hall  (1988).  However,  in  situations  in  which  the  Gortler  number  is  large, 
the  results  of  this  work  would  be  of  limited  relevance  since  we  would  expect  that  in  this 
case  the  most  unstable  linear  Gortler  mode  would  soon  dominate  the  flow  field.  In  fact, 
the  receptivity  problem  for  the  most  unstable  mode  has  been  considered  by  Denier  et  al 
(1991)  who  demonstrated  that  the  coupling  coefficient  between  a  surface  perturbation  and 
the  induced  velocity  field  is  an  0(1)  quantity  and  as  such  we  would  expect  these  modes  to 
be  generated  by,  for  example,  isolated  roughness  elements  on  the  surface. 

To  consider  the  nonlinear  development  of  the  most  unstable  Gortler  mode  we  define  new 
variables 

r  =  G'/"y,  Z  =  G’/^z,  X  =  G"/"(x-x*),  (2.6) 

and  expand  the  velocity  field  and  pressure  in  the  form 

(u,  V,  W,p)  =  G*/"(G-"/"t/,  F,  W,  +  ■  •  • .  (2.7) 

Here,  the  scaling  in  X  anticipates  the  result  from  Denier  et  al  (1991)  that  the  downstream 
growth  rate  of  the  most  unstable  Gortler  mode  is  0(G^/^).  The  governing  equations  (2.4) 
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then  become,  to  leading  order  in  powers  of  G 


Gx  +  Vy  +  Wz  =  0, 

Vyy  +  Vzz  -  =  ixYVx  +  UUx  +  VUy  +  WUz. 

VVy  +  -  x{x*)lJyU  -Py  =  hYVx  +  UVx  +  VVV  +  WVz  + 

Wyy  +  Wzz  -Pz  =  f^YWx  +  UWx  +  VWy  +  WWz, 

where  we  have  made  use  of  the  fact  that  u  ~  fi(x*)y  +  •  •  •  as  y  — »  0.  With  a  suitable 
rescaling  we  can  remove  the  constant  coefficients  in  the  above  system,  and  hence,  without 
loss  of  generality,  we  will  assume  ^  =  x  =  1  in  the  above  system.  The  governing  equations 
are  then 

Ux  +  Vy  +  Wz  =  0, 

Uyy  +Uzz-V  =  YUx  +  UUx  +  VUy  +  WUz, 

1  2  (2.8) 

VVy  +  VVz  -YU -Py  =  YVx  +  UVx  +  VVV  +  WVz  +  ^U\ 

Wyy  +  Wzz  -Pz  =  YWx  +  UWx  +  VWy  +  WWz. 

The  boundary  conditions  appropriate  to  the  above  system  are  found  to  be 

U=-V  =  W  =  0  F  =  0, 

(2.9) 

U  A{X)  r  ->  oo, 


where  the  displacement  function  A  is  independent  of  Z  (ensuring  that  the  nonlinear  vortex 
state  is  confined  to  the  wall  layer).  To  complete  the  description  of  the  flow  regime  wc 
require  a  region  where  y  =  0(1)  in  which  the  mean  flow  adjusts  to  its  free  stream  value. 
However  that  region  is  passive  so  we  shall  not  consider  it  here. 

In  order  to  reduce  the  system  (2.8)  to  a  form  suitable  for  computational  purposes  we 
eliminate  W  and  P  from  the  linear  terms  in  (2.8c, d)  to  give 

^  1^)  (2.10) 

where 

iVi  =  UUx  +  VUy  +  WUz, 

=  UVx  +  H  W  +  WVz  +  lu\  (2.11) 

Ns  =  UWx  +  VWy  +  WWz. 
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The  numerical  procedure  employed  in  the  solution  of  (2.8a,b)  and  (2.10)  subject  to  the 
boundary  conditions  (2.9)  is  based  on  that  used  by  Hall  (1988),  in  the  solution  of  the 
nonlinear  Gdrtler  equations  at  0(1)  Gdrtler  numbers,  and  for  this  reason  we  give  only  the 
salient  details  of  the  scheme  here  (the  reader  is  referred  to  the  aforementioned  paper  for 
full  details). 

.Anticipating  the  well  known  result  that  nonlinear  interactions  in  the  Gdrtler  problem  do 
not  generate  a  mean  flow  component  in  the  spanwise  direction  we  write 

OO  OO 

{U,V)  =  {Uo,Vo)  +  Vn)  cos  nkz,  W  =  '^Wnsinnkz,  (2-12) 

n=l  Ti=l 

where,  without  loss  of  generality,  we  have  chosen  the  origin  in  the  spanwise  direction  such 
that  U  and  V  are  even  in  Z  whilst  W  is  odd  in  Z.  Substituting  the  expansions  (2.12) 
into  (2.8a,b)  and  (2.10)  and  equating  like  Fourier  components  gives,  for  the  mean  flow 
correction 

Uqyy  ■■  1^0  —  y^Uox  4*  UoUox  "b  VqUqy  "b  Fo,  (2.13) 


where 


1  ^ 

=  9  H  “  UmVmY  -  2kmUmWm)  , 


whilst  Vo  is  determined  from 


dUo  dVo_ 

dX  dY 


(2.14) 


From  a  computational  standpoint  we  must  necessarily  truncate  the  infinite  sums  occurring 
in  (2.12)  at  some  suitably  large  value  of  the  upper  limit.  Therefore  we  will  formally  replace 
the  upper  limit  in  (2.12)  by  N. 

The  governing  equation  for  the  higher  harmonics  terms  is  then  given  by 

Ur^yy  -  U^k^Un  "  (1  +  UoyWu  -  {Y  +  Uo)U„X  =  Fn 

N-\ 

—  ^  ^  Yn  —  mUmY  Un  —  mYmY  —  ^^Un  —  mY^rn 


"b  ^  ^  Yn-\-mFmY  Un-^-rnYmY  ^^YWn-^rnUm  n-\- m^Y n 


(2.15) 


N 

"b  ^  ^  ^m  —  n^'niY  F  m  —  nYmY  — nFm  ^kUm  —  n^Ym  i 

m  =n  +  I 
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where  we  have  made  use  of  the  continuity  equation  to  eliminate  the  X-derivatives  on 
from  the  nonlinear  terms  on  the  right  hand  side  cf  (2.15).  A  similar  equation  governing 
the  evolution  of  the  components  Vn  can  be  obtained  from  (2.10)  which,  for  the  sake  of 
brevity,  is  not  presented  here. 

We  now  proceed  with  a  brief  description  of  the  scheme  used  to  solve  (2.13),  (2.15)  together 
with  the  corresponding  equation  for  the  V„.  Suppose  then  that  Uq,  Vq  and  Z7„,  F„,  W„,  for 
n  =  1,  2,  •  ■  •  ,  X  are  known  at  some  particular  X-station,  Xq  say.  The  mean  flow  equation 
(2.13)  is  discretized  using  finite  differences  in  the  X-  and  K-directions  according  to 


Here  the  indices  m  and  n  refer  to  the  grid  point  X  =  Xq  +  m€,Y  =  nh.  Initially  the 
nonlinear  terms  on  the  right-hand  side  of  (2.16)  are  evaluated  with  k  =  m-l  which  yields  a 
tridiagonal  system  for  Uo  a.t  X  =  Xo  +  (m4-l)e.  Similarly  the  finite  difference  discretization 
of  (2.15)  can  be  used  to  give  Um  (m  =  1,---,X)  at  X  =  Xq  +  (m  +  l)e  as  can  the 
corresponding  equation  for  Vm,  but  now  by  solving  a  pentadiagonal  system.  The  nonlinear 
terms  can  now  be  expressed  in  terms  of  the  velocity  field  evaluated  at  X  =  Xo  +  (m  +  l)e 
and  the  equations  are  again  solved  for  the  flow  quantities  at  X  =  Xq  +  (m  +  1  ''e,  repeating 
the  iteration  procedure  until  the  changes  in  etc  are  sufficiently  small.  In  this  way 

(2.16)  and  the  corresponding  equations  for  Um-,ym  are  solved  for  fc  =  m  by  iterating  on 
the  nonlinear  terms  on  the  right-hand  sides.  This  scheme  can  then  be  used  to  march 
the  solution  downstream  from  some  initial  position  given  some  suitable  form  of  initial 
perturbation  to  the  velocity  field. 

3.  Results  and  discussion 

We  now  describe  the  results  obtained  from  the  numerical  scheme  outlined  above.  All 
results  presented  are  for  the  choice  of  wavenumber  k  =  0.476  which  corresponds  to  the 
critical  wavenumber  of  the  most  unstable  Gdrtler  mode  identified  by  Denier  et  al  (1991); 
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some  discussion  on  the  form  of  the  solutions  for  different  values  of  the  wavenumber  will  be 
given.  The  basic  state  was  perturbed  at  X  =  1  by  imposing  the  disturbance 

exp(-yV2),  F,(r)  =  o,  (3.1) 

and  integrating  the  linearized  equations  to  X  =  101.  Here  A  is  a  measure  of  the  disturbance 
amplitude.  At  this  stage  the  disturbance  has  attained  its  maximum  growth  rate  according 
to  linear  theory  (see  Denier  et  al  (1991)).  The  linear  velocity  field  was  then  given  an 
amplitude  equal  to  the  maximum  X-velocity  component  whilst  all  higher  harmonics  were 
set  equal  to  zero.  The  nonlinear  equations  were  then  integrated  for  X  >  101. 

We  note  that  the  precise  form  of  the  initial  perturbation  (3.1)  is  of  little  importance  in  the 
linear  stage  of  the  evolution  of  the  most  unstable  mode.  This  can  be  seen  by  considering 
the  linearized  version  of  the  governing  equations;  a  comparison  between  the  present  work 
and  the  eigenvalue  problem  of  Denier  et  al  (1991)  (see  their  equation  (5.14))  shows  that  the 
linearized  governing  equations,  subject  to  the  initial  condition  (3.1),  constitute  an  initial 
value  problem  corresponding  to  the  given  spanwise  wavenumber.  Thus,  given  any  reason¬ 
able  initial  condition  (ie  one  which  satisfies  ?U  the  boundary  conditions  on  the  disturbance 
velocity  field)  the  solution  of  the  linearized  equations  will  evolve  downstream  into  the  char¬ 
acteristic  shape  of  the  eigenfunction  of  Denier  et  al  (1991)  (in  fact,  the  downstream  form 
of  the  velocity  field  will  be  a  constant  multiple  of  this  eigenfunction)  and  the  linear  growth 
rate  will  approach  the  value  found  from  the  eigenvalue  problem;  for  the  results  presented 
here  this  corresponds  to  k  =  0.476,  /?  =  0.312  (here  (3  is  the  linear  growth  rate).  This 
should  be  compared  with  the  corresponding  results  of  Hall  (1983),  for  the  evolution  of 
vortices  with  0(1)  Gortler  number,  where  the  growth  rate  and  the  structure  of  the  linear 
modes  exhibit  a  sensitive  dependence  on  both  the  form  of  the  initial  disturbance  and  the 
position  at  which  the  initial  disturbance  is  imposed.  This  difference  is  readily  explained  by 
the  fact  that  the  most  unstable  mode  of  Denier  et  al  (1991)  evolves  in,  what  is  essentially,  a 
parallel  manner  (due  in  turn  to  the  small  x  scale  over  which  these  modes  develop)  whereas 
the  linear  modes  of  Hall  (1983  )  exhibit  an  inherent  nonparallel  structure  (due  to  tlie  0(1) 
j-  scale  over  which  these  modes  develop). 
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In  order  to  measure  the  growth  of  the  harmonic  components  of  the  flow  field  we  define  the 
energy  of  the  nth  harmonic  to  be 

OO 

En  =  J  {U^iX,Y)  +  V^iX,Y)  +  W^{X,Y)}dY,  n  =  (3.2) 

0 

Note  that  we  do  not  calculate  the  energy  of  the  mean  flow  correction  as  in  Hall  (ly88) 
since  in  the  present  case  we  have  Uq  — >  Uq{X)  v  hen  Y  — >  oo.  The  growth  rate  <7„(X)  of 
the  nth  mode  is  then  defined  as 


In  the  absence  of  the  nonlinear  terms  the  growth  rate  (X)  would  be  twice  that  calculated 
by  Denier  et  al  (1991);  in  light  of  the  discussion  of  the  previous  paragraph  (7i{X)  0.624 

as  A’  ^  oo  (see  figure  2). 

With  the  perturbation  to  the  basic  state  given  by  (3.1)  the  governing  equations  were  solved 
f<>r  the  choice  of  parameters  h  —  0.1,  e  =  0.001  with  Foo  =  100  (corresponding  to  1000  grid 
points  in  the  normal  direction).  These  values  were  chosen  after  significant  testing  of  the 
code  for  accuracy  ai  d  stability.  With  h  fixed  the  calculations  were  carried  out  for  various 
values  of  e  until  the  results  from  consecutive  runs  agreed  to  within  graphical  accuracy.  The 
otreamwise  step  size  c  was  then  fixed  and  h  was  decreased,  with  Yoo  fixed,  until,  again,  the 
results  from  consecutive  calculations  had  converged,  to  within  graphical  accuracy.  We  then 
fixed  both  e  and  h  and  varied  Y^o  which  resulted  in  little  or  no  change  in  the  computed 
flow  quantities.  It  was  found  that  the  numerical  scheme  was  stable  for  e  =  O(h^).  In  figure 
3  we  present  a  comparison  between  calculations  performed  for  e  =  0.002  and  e  =  0.001 
with  the  normal  step-length  /i  =  0.1  fixed  (see  figure  caption).  Up  until  the  stage  at 
which  the  calculation  for  e  =  0.001  fails  the  results  are  identical  (note  that  the  results  for 
e  =  0.002  have  been  offset).  However,  there  is  a  small,  yet  significant,  discrepancy  on  th“ 
position  at  which  the  breakdown  occurs.  We  interpret  this  discrepancy  in  the  breakdown 
position  which  occurs  for  the  larger  of  the  two  values  of  e  as  being  due  to  the  numerical 
scheme,  essentially,  jumping’  over  the  breakdown  position  and  continuing  into  a  region 
where,  in  fact,  the  numerical  code  is  no  longer  valid.  Various  calculations  were  done  for 
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which  =  4,8  and  16  Fourier  modes  in  thr*  spanwise  direction  were  retained;  the  results 
presented  are  for  the  case  N  =  8  (the  re-  Its  from  the  calculations  for  =  4  and  A^  =  16 
show  little,  if  any  variations  from  those  obtained  for  N  =  8).  In  figure  4  we  present  a 
comparison  of  the  results  obtained  with  AT  =  8  and  A^  =  16  Fourier  modes  retained  in 
the  expansions  (2.12).  Again  the  results  are  identical,  to  within  graphical  accuracy  foi 
the  two  cases.  (Note  that  the  discrepancy  in  the  breakdown  positions  for  the  N  =  8  and 
A^  =  16  calculations  is  due  to  the  increased  streamwise  step-length  required  to  make  the 
latter  calculation  feasible).  The  parameter  values  quoted  above  were  then  chosen  with 
two  criteria  in  mind;  firstly  they  were  chosen  to  be  suitably  small  in  order  to  resolve  both 
the  normal  and  streamwise  structure  of  the  flow  in  the  later  stages  of  the  calculations 
and  secondly  were  chosen  so  as  to  allow  the  calculations  to  be  performed  in  an  acceptable 
amount  of  time.  All  results  presented  were  obtained  on  a  Cray  Y-MP  situated  at  the 
NASA  Langley  Research  Centre. 

In  figure  5  we  present  a  plot  of  the  energy  of  the  nth  harmonic  as  a  function  of  downstream 
position  X.  Here  we  see  that,  after  an  initial  interval  of  decay  due  to  the  solution  rapidly 
adjusting  from  its  linear  state  to  a  fully  nonlinear  state  in  which  the  mean  flow  correction 
is  driven  by  the  vortex  velocity  field,  the  energy  in  each  mode  is  a  monotonically  increasing 
function  of  x.  This  was  found  to  continue  until  such  point  as  the  energy  in  each  harmonic 
undergoes  a  period  of  rapid  growth  (occurring  over  a  small  streamwise  interval).  However, 
at  this  stage  of  the  calculation  a  jjoint  of  separation  has  been  encountered  resulting  in 
reversed  flow.  Our  numerical  scheme  is  then  no  longer  valid,  as  it  relies  on  the  parabolic 
nature  of  the  governing  equations  and  as  such  is  not  able  to  deal  with  the  reversed  flow 
in  the  streamwise  direction.  Numerous  calculations  were  performed  in  order  to  verify  this; 
in  particular,  increasing  the  number  of  modes  retained  in  our  expansions  from  N  —  8  to 
N  -  16.  The  results  from  this  calculation  were  found  to  agree  with  that  for  Y  =  8  Fourier 
modes  to  within  the  graphical  accuracy  of  figure  5.  Similarly,  decreasing  the  step  size  in 
the  streamwise  direction  was  found  to  have  little  effect  on  the  results  of  the  calculation. 
All  results  presented  were  obtained  for  the  input  disturbance  amplitude  A  =  0.075.  In 
order  to  test  the  effect  of  varying  the  amplitude  of  the  initial  disturbance  velocity  field 
we  repeated  our  calculations  fo’  various  values  of  the  amplitude  A.  The  results  of  these 
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calculations  demonstrated  that  varying  A  had  no  qualitative  effect  upon  the  breakdown 
which  is  encountered  for  the  case  presented,  namely  A  =  0.075;  decreasing  A  has  the 
effect  of  delaying  the  breakdown.  For  this  reason  we  restrict  our  attention  to  the  case 
A  =  0.075;  the  results  presented  can  then  be  considered  as  typical  for  the  whole  range  of 
input  amplitudes. 

Before  proceeding  with  a  discussion  of  this  breakdown  we  make  some  comparison  of  our 
results  with  those  of  Hall  (1988).  In  figure  6  we  present  a  plot  of  the  growth  rates  cr„,  (n  = 
1,  •  •  •  ,8)  as  a  function  of  the  downstream  position  x.  Here  we  see  that  the  growth  rates, 
after  an  initial  adjustment  stage  wherein  the  flow  adjusts  in  order  to  accommodate  the 
nonlinear  terms,  increase  with  increasing  x;  in  fact  the  growth  rates  are  ordered  such  that 
tjj  >  CTj-i,  {j  =  2,  •  •  • ,  8).  However,  the  growth  rate  of  the  first  harmonic,  a\ ,  is  decreased 
from  that  found  in  the  absence  of  the  nonlinear  terms;  in  this  respect  nonlinear  effects 
are  stabilizing.  This  should  be  compared  with  the  results  of  Hall  (1988)  who  found  a 
similar  result  in  that  the  growth  of  the  first  harmonic  is  decreased  in  comparison  with  that 
obtained  from  the  linearized  system.  However,  an  important  distinction  arises  between 
the  present  work  and  that  of  Hall  (1988)  in  that  for  the  case  of  0(1)  Gortler  numbers  and 
wavenumbers  the  energy  is  distributed  predominantly  in  the  first  harmonic  and  mean  flow 
correction  whereas  in  the  present  case  no  such  partition  of  energy  is  found. 

We  now  turn  our  attention  to  a  discussion  of  the  breakdown  alluded  to  above.  As  noted 
earlier  our  numerical  scheme  was  found  to  breakdown  at  some  downstream  position.  This 
breakdown  was  characterized  by  the  normal,  and  subsequently  the  spanwise,  components 
of  the  velocity  field  undergoing  a  period  of  rapid  growth  over  a  small  streamwise  inter¬ 
val,  typically  of  the  order  of  5  streamwise  steplengths,  whereas  the  streamwise  velocity 
component  shows  little,  if  any,  variation  over  this  interval. 

In  order  to  determine  whether  this  breakdown  was  due  to  a  singularity  of  the  governing 
equations,  we  repeated  our  calculations  with  an  increased  number  of  Fourier  modes  re¬ 
tained.  As  each  higher  harmonic  is  forced  by  its  lower  harmonics  we  would  expect  that, 
if  such  a  singularity  were  encountered,  increasing  the  number  of  Fourier  modes  retained 
in  our  calculations  would  have  a  dramatic  effect  upon  the  latter  stages  of  the  calculation. 
However,  this  wa»s  found  not  to  be  the  case.  The  inclusion  of  additional  Fourier  modes  was 
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found  to  have  little  effect  on  the  breakdown.  For  this  reason  we  feel  justified  in  eliminating 
the  possibility  of  a  spanwise  localized  singularity  in  the  governing  equations  as  being  the 
underlying  cause  for  the  breakdown  which  is  encountered. 

In  light  of  the  above  discussion  an  alternative  mechanism  for  the  breakdown  must  be  found. 
An  obvious  candidate  immediately  presents  itself  in  the  form  of  a  singularity  due  to  flow 
separation.  From  the  work  of  Hall  (1988)  we  know  that  flow  separation  is  encountered  in 
the  problem  of  0(1)  wavelength  Gortler  vortices  at  0(1)  Gortlcr  numbers  (see  also  the 
discussion  in  Hall  and  Horseman  (1990)).  As  such,  flow  separation  would  appear  as  the 
likely  candidate  for  the  breakdown  experienced  in  the  present  problem. 

In  order  to  test  this  conjecture  we  calculated  the  total  skin  friction 

N 

tt(X,Z)  =  1  +  ^  (X,0)  cos  A:nZ,  (3.4) 

n=0 

at  each  downstream  position  on  a  uniform  spanwise  grid  of  100  points  over  one  half  vortex 
wavelength.  Here,  we  note  that  the  first  term  on  the  right  hand  side  of  (3.4)  arises  from 
the  uniform  shear,  u  ~  Y,  in  the  wall  layer  (see  equation  (2.S)).  The  point  (X*,Z*)  of 
flow  separation  is  then  defined  by 


rT(X*,Z*)^0. 

At  such  a  point  our  finite  difference  discretization  in  the  downstream  variable  X  is  no 
longer  valid  as  it  does  not  allow  for  reversed  flow.  In  figure  7  we  present  a  contour  plot  of 
(tt  -  1)  in  the  (AT,  Z)-plane  for  X  in  the  range  (108, 108.037)  and  Z  in  the  range  (0,  tt/o). 
Here  we  see  that,  after  an  interval  of  little  variation,  the  total  skin  friction  rapidly  decreases 
until  the  last  downstream  location  at  which  (tt—1)  is  negative  and  separation  has  occurred. 
At  this  point  in  our  calculation  the  numerical  scheme  used  is  no  longer  valid  for  the  reason 
given  above.  We  found  it  impossible  to  accurately  resolve  the  position  at  which  separation 
occurs,  though  not  through  any  deficiency  in  the  code,  but  due  to  the  fact  that  we  had 
reached  the  limit  of  our  available  computing  resources. 
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Not  surprisingly,  the  point  of  separation  is  found  to  occur  at  the  position  Z  =  tt/A;,  the 
position  at  which  upwelling  occurs.  In  figure  8  we  present  a  series  of  contour  plots  of 
constant  total  streamwise  velocity 

N 

Ut{X,  Y,Z)=Y  +  Yi  Unix,  Y)  cos  knZ, 

n=0 

at  various  downstream  positions  which  shows  the  Initial  growth  of  the  vortex.  In  the 
absence  of  the  nonlinear  vortex  these  plots  vvould  appear  as  a  series  of  equally  spaced 
horizontal  lines.  The  growth  of  the  vortex  is  then  represented  by  increased  curvature  of 
the  streamlines  of  constant  streamwise  velocity.  At  X  —  102.0  (figure  8a)  the  nonlinear 
vortex  has  begun  to  grow  until  at  A  =  108.0  (figure  8d)  we  see  the  characteristic  shape 
for  Gortler  vortices  (see  Hall  and  Horseman  (1990)). 

In  figure  9a-e  we  again  present  contour  plots  of  the  total  streamwise  velocity  at  the  last 
five  streamwise  locations  preceding  the  breakdown  of  our  numerical  code.  At  this  stage 
we  see  the  development  of  a  secondary  structure  localized  about  Z  =  n /a  m  the  vicinity 
of  the  wall  (see  the  highlighted  region  in  figure  9a).  This  structure  increasingly  distorts 
the  underlying  streamlines  until  at  the  last  point  of  our  calculation  (figure  9e)  the  flow 
is  about  to  become  reversed.  As  noted  above  we  were  unable  to  continue  our  calculation 
past  the  point  given  for  figure  9e,  however  we  would  expect  that  if  such  a  calculation  were 
possible  we  would  find  a  small  region  of  reversed  flow  centered  about  the  point  Z  —  -k  j a 
which  would  subsequently  increase  in  size  as  the  calculation  were  continued  downstream. 
4.  Conclusion 

We  have  investigated  the  spatial  nonlinear  evolution  of  the  fastest  growing  Gortler  vortex 
which  can  occur  in  a  boundary  layer  flow.  The  mode  occurs  in  a  quasi-parallel  form  only 
at  high  Gortler  numbers  so  our  investigation  is  particularly  relevant  to  flows  where  high 
Reynolds  numbers  and  moderate  curvatures  are  possible.  Hence  our  investigation  is,  for 
example,  relevant  to  the  flow  over  turbine  blades  and  also  to  the  flow  in  jet  engine  inlets 
though  compressible  effects  would  certainly  be  present  in  the  latter  case.  A  crucial  feature 
of  the  mode  is  that  it  occurs  adjacent  to  the  wall  so  that  relatively  high  shear  stresses  will 
be  induced  by  its  presence. 
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Our  calculations  suggest  that  the  nonlinear  development  of  the  most  unstable  Gortler 
mode  results  in  a  reversed  flow.  Although  such  a  point  of  separation  gives  little  hope 
of  following  the  development  of  these  modes  downstream  past  the  point  of  separation 
we  note  one  important  point  which  arises  from  our  calculations.  Recent  work  on  the 
secondary  instability  of  nonlinear  Gortler  vortices  at  0(1)  Gortler  number  by  Hall  and 
Horseman  (1990)  has  demonstrated  that  the  fully  nonlinear  vortex  velocity  fields  obtained 
numerically  by  Hall  (1988)  are  highly  susceptible  to  an  inviscid  Rayleigh  instability.  This 
mode  of  instability  is  particularly  relevant  to  flows  for  which  the  streamwise  velocity  field 
becomes  inflexional  in  both  the  streamwise  and  spanwise  directions.  Such  a  situation  arises 
in  the  present  problem  and  hence  we  would  expect  the  flow  to  be  radically  altered,  due 
to  the  linear  growth  of  Rayleigh  modes  and  their  subsequent  nonlinear  evolution,  before 
the  point  of  separation  is  approached,  as  is  the  case  for  Gortler  vortices  at  0(1)  Gortler 
numbers  (see  Hall  and  Horseman  (1990)  for  a  discussion)).  This  question  is  currently 
under  investigation. 

Thus  we  conclude  then  by  noting  that  our  calculations  suggest  that  in  relatively  highly 
curved  boundary  layers,  transition  will  be  induced  by  the  Rayleigh  instability  breakdown  of 
the  vortex  structures  driven  by  the  Gortler  instability.  It  should  be  noted  that  alternative 
instabilities  such  as  Tollmien-Schlichting  waves  do  not  have  growth  rates  as  large  as  those 
of  Gortler  vortices  or  Rayleigh  waves  and  are  therefore  likely  to  be  unimportant  in  the 
transition  process. 
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Figure  1.  A  schematic  representation  of  the  growth  rate  versus  wavenumber  in  the  large 
Gortler  number  limit.  Marked  are  the  regions  corresponding  to  the  inviscid,  most  unstable 
and  right-hand  branch  Gortler  modes,  respectively. 


Figure  4.  Comparison  of  results  obtained  by  retaining  different  number  of  Fourier  modes. 
Solid  line  iV  =  8,  dashed  line  N  =  16.  In  order  to  distinguish  the  plots  for  =  16  these 
have  been  offset  by  —0.5.  Due  to  limitations  on  available  computing  time  the  results  for 
N  =  16  were  obtained  for  a  streamwise  step-length  of  e  =  0.01. 
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Figure  9.  Contour  plot  of  the  total  downstream  velocity  field  in  the  final  stages  of  the 
calculation;  (a)  X  =  108.032  (b)  X  =  108.033,  (c)  X  =  108.034,  (d)  X  =  108.035,  (e) 
X  =  108.036. 
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The  nonlinear  development  of  the  most  unstable  Gortlervortex  mode  in  boundary  layer  flows  over  curved  walls 
is  investigated.  The  most  unstable  Gortlermode  is  confined  to  a  viscous  wall  layer  of  thickness  and  has 

spanwise  wavelength  0(G^*7*);  it  is,  of  course,  most  relevant  to  flow  situations  where  the  Gortlernumber  G  »  1. 
The  nonlinear  equations  governing  the  evolution  of  this  mode  over  an  0(G‘~^f^)  streamwise  lengthscale  are  derived 
and  are  found  to  be  of  a  fully  nonparallel  nature.  The  solution  of  these  equations  is  achieved  by  making  use  of 
the  numerical  scheme  used  by  Hall  (1988)  for  the  numerical  solution  of  the  nonlinear  Gortlerequations  valid  for 
0(1)  Qortlernumbers.  Thus,  the  spanwise  dependence  of  the  flow  is  described  by  a  Fourier  expansion  whereas  the 
streamwise  and  normrd  variations  of  the  flow  are  dealt  with  by  employing  a  suitable  finite  difference  discretization 
of  the  governing  equations.  Our  calculations  demonstrate  that,  given  a  suitable  initial  disturbance,  after  a  brief 
interval  of  decay,  the  energy  in  all  the  higher  harmonics  grows  until  a  singularity  is  encountered  at  some  downstream 
position.  The  structure  of  the  flow  field  as  this  singularity  is  approached  suggests  that  the  singularity  is  responsible 
for  the  vortices,  which  are  initially  confined  to  the  thin  viscous  wall  layer,  moving  away  from  the  wall  and  into  the 
core  of  the  boundary  layer. 
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